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DIFFUSION  APPROXIMATIONS  FOR  THE  COOPERATIVE  SERVICE 
OF  VOICE  AND  DATA  MESSAGES 


by 


J.  P.  Lehoczky 
Carnegie-Mellon  University 
Pittsburgh,  PA 

and 

D.  P.  Gaver 
Naval  Postgraduate  School 
Monterey,  CA 


INTRODUCTION 

In  this  paper  we  study  the  behavior  of  a  queueing 
system  which  arises  in  the  study  of  certain  communication 
networks.   Specifically  we  study  a  queueing  phenomenon  which 
arises  with  the  SENET  network,  as  described  by  Coviello 
and  Vena  (1975)  or  Barbacci  and  Oakley  (1976)  .   This  network 
allows  for  both  voice  and  data  messages  to  be  transmitted 
over  the  same  channels  by  using  a  special  type  of  integrated 
circuit  and  packet-switched  multiplexor  structure.   The  two 
classes  of  traffic  have  substantially  different  performance 
requirements.   Voice  messages  tend  to  possess  great  redundancy, 
and  hence  not  to  be  sensitive  to  channel  error  rates,  while 
data  is  very  sensitive  to  channel  error,  having  essentially 
no  redundancy.   Voice  messages  on  the  other  hand  have  critical 
timing  requirements  and  cannot  be  queued,  while  data  is 


relatively  insensitive  to  timing  and  can  be  queued.   Addition- 
ally, voice  messages  tend  to  be  very  long  relative  to  data 
messages  which  can  be  broken  up  into  small  packets.   These 
special  requirements  have  led  to  the  following  queueing  net- 
work.  A  node  of  the  network  consists  of  c  +  v  channels  or 
servers.   The  voice  messages  are  assigned  to   v   channels  and 
do   not  queue.   Thus  the  voice  messages  operate  as  a  loss 
system.   Data  messages  may  use   c   channels  exclusively  and 
any  unused  voice  channels;  however,  voice  preempts  data  using 
voice  channels.   Data  messages  are  queued  if  necessary. 
Typical  performance  measures  that  one  may  wish  to  calculate 
include  the  loss  rate  of  voice  traffic  and  the  mean  data 
queue  length. 

We  make  standard  probabilistic  assumptions.   Specifically, 
we  assume  voice  traffic  arrives  according  to  a  Poisson(A) 
process  and  each  voice  message  has  an  independent  exponential ( p) 
service  time.   Data  messages  are  assumed  to  have  independent 
exponential (n)  service  times  and  arrive  according  to  a  Poisson(6) 
process.   With  these  assumptions  voice  is  an  M/M/v/v  loss 
system,  and  data  is  an  M/M/S  system  where   S=c+v-V(t) 
with   V(t)  =  number  of  voice  messages  in  service.   The 
stochastic  process   { (X (t) , V( t) ) ,  t  >  0}   is  Markov  with 
state  space   Z   x  {0,1,..., v}   where   X(t)  =  data  system  size 
at  time  t.   One  can  easily  write  the  Kolmogorov  forward 
equations  appropriate  for  this  system;  however,  these  equations 
do  not  yield  a  closed  form  solution.   To  describe  this  system 
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one  must  either  numerically  solve  the  forward  equations  or 
introduce  approximations . 

This  system  has  been  studied  previously  by  a  number 
of  researchers  including  Halfin  and  Segal  (1972) ,  Halfin  (1972) , 
Fischer  and  Harris  (1976),  Bhat  and  Fischer  (1976),  Fischer 
(1977) ,  Chang  (1977) ,  and  Gaver  and  Lehoczky  (1979a, b) .   The 
last  two  papers  introduce  a  "fluid  flow"  and  a  diffusion 
approximation  and  derive  explicit  formulas  for  data  queue 
behavior.   These  papers  focus  on  the  important  case  in  which 
p ,  =  6 /n  >  c.   In  such  a  situation  the  data  messages  must 

have  access  to  voice  channels  for  the  system  to  be  stable. 

4 
Furthermore,  it  was  assumed  that  r\/\i      was  large,  say  10  . 

Under  these  circumstances  the  data  flow  could  be  treated 

deterministically .    Suppose  we  define   p   =  A/y   and 

q  =  (Pv/v! )/£._«  P  / j ! ,    the  Erlang  B  blocking  probability. 

The  total  traffic  intensity  on  the  c  +  v  channels  is  given 

by   p,  +  p  (1-q) ,  or  we  could  define   p  =  (p,  +  p  (1-q) ) / (c+v) . 

A  heavy  traffic  approximation  can  be  derived  for  this  case 

p  f     1.   Such  an  approximation  was  derived  in  Gaver  and 

Lehoczky  (19  79b)  assuming  x\/\x      was  large;  a  Wiener  process 

with  reflecting  boundary  was  found  appropriate.    In  this 

paper  we  derive  a  heavy  traffic  approximation  for  the  system 

without  the  fluid  flow  assumption  that   n/y   is  large. 

The  methodology  is  drawn  heavily  from  the  approach  of  Burman 

(1979) .   In  this  approach  one  characterizes  a  Markov  process 


by  its  infinitesimal  generator.   One  next  suitably  normalizes 
the  process  so  that  the  generator  converges  to  a  limiting 
infinitesimal  generator  (in  this  case  to  that  of  a  reflected 
Brownian  motion) .   This  convergence  allows  the  conclusion  that 
the  finite  dimensional  distributions  of  the  normalized  Markov 
process  converge.    The  diffusion  approximation  consists  of 
treating  the  actual  process  through  its  limiting  behavior. 
The  details  are  somewhat  complicated  by  the  presence  of  a 
boundary . 


2. 

Let   {(X(t),V(t)),  t  _>  0}   be  a  bivariate  Markov 
process  with  state  space   S  =  Z   x  {0,1,... ,  v)  •   Here 
{V(t),  t  _>  0}   is  marginally  an  M/M/v/v  loss  system  with 
arrival  rate   A   and  service  rate   y.   Conditional  on   V(t) , 
{X(t),  t  >_  0}   is  an  M/M/(c  +  v  -  V(t))   queueing  system  with 
arrival  rate  6   and  service  rate   n .   We  say  that  the   V 
process  subordinates  the   X   process.   We  let 


s  = 


tt 


v 

-d+pv) 


V 

-<2+p  ) 


V 


v-1 


V 

(v-l+p  ) 


(2.1) 


the  infinitesimal  generator  of  the   V  process. 

The  generator  of  the   (X,V)   process  is  given  by 


Af(x,k)  = 


Qf(x,k)  +  <5(f(x+l,k)  -  f(x,k)) 

+  n(c+v-k) (f (x-l,k)-f (x,k)) 
if   x  2  c+v-k 

Of (x,k)  + 6  (f (x+l,k)  -  f(x,k)) 

+  nx(f  (x-l,k)  -  f  (x,k)  ) 

if   x  <  c+v-k 


(2.2) 


for   f:S  ■*■  R  continuous  where 


Qf(x,k)  =  p  f(x,k+l)  -  (k+pv)  f  (x,k)  +  kf  (x,k-l) 


v  >  k  >  0 


(2.3) 


and   f(x,-l)  =  f(x,v+l)  =  0.   Clearly   Qf(x)  =  0,  that  is   Q 
annihilates  functions  of   x   alone.   We  next  normalize  the 
(X,V)  process  by  defining   X  (t)  =  X(nt)//n   and  V  (t)  =  V(nt) 
One  can  calculate  the  generator   of  the  Markov  process 
{ (X  (t) ,  V  (t) ) ,  t  >  0}   having  state  space 


n 


n 


S   =  {0,  l//n,  2//n,  .  .  .  }  x  {0,1,  .  .  .  ,v}   to  be 


/ 


nQf(x,k)  +6n(f(x+  l//~h,k)  -  f(x,k)) 


+  nn(f  (x  -  1/Tn,k) 


C  +  V  - 


if   x  > 


A  f  (x  k)  =J 

nQf(x,k)  +5n(f(x  +  l//n,k) 


f (x,k)) 
k 


-  f (x,k) ) 


(2.4) 


+  nn  /nx  (f(x  +  l//n,k)  -  f(x,k)) 


if   x  =  0,  l//n,...,  (c+v-k)//n" 


We  assume  f(x,k)  has  three  bounded  derivatives  in  x 
for  each  fixed  k.  With  this  assumption  one  can  expand  terms 
in  (2.4)  in  a  Taylor  series  and  rewrite  as 
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r 


nQf(x,k)    +  n1/2f    (x,k) (6 -n (c+v-k) ) 


+    2    fxx(x'k)  (5+rl(c+y-k)  ) 


+   Otn"1/2) 


if      x    >    (c+v-k) //h 


A  f(x,k) 


=< 


nQf(x,k)    +   n1/2f    (x,k)  (6  -n/nx) 


+   \  fxx(x,k)  (5+n   /nx)    +   0(n   1/2) 


if      x   =    0,    l//n, 


,     (c+v-k) //n 


(2.5) 


with      f    (x,k)    =  ~  f  (x,k)      and      f      (x,k)    =   -^f(x,k). 

a  X 

We    ultimately  wish    to   prove    that    the    finite    dimensional 

distributions   of      (X    (t)  ,    t  ^  0}      converge    to    those   of   a 

Wiener   process   with   reflecting  barrier   at   the   origin.      This 

can  be    restated   in    terms    of    semi-groups.      We    let      {Tj_/    t  ^  0} 

be  the  semi-group  of  operators  associated  with 

{ (X  (t) ,  V  (t)),  t  >  0}   and   {T\,  t  >  0}   be  that  associated 
n      n        —  t    — 

with  a  Wiener  process  having  reflecting  barrier  at  0.   Let 
g  be  a  continuous  function   g:R'  ■*■   R'  .   Knowledge  of  the 
semi-group  is  equivalent  to  knowledge  of  the  transition 

functions  by  taking  a  sequence  of   g's  which  approximate 

i  n  oo      | 

indicator  functions.   We  wish  to  prove   I T   g(x,k)  -  T   g(x)|  -*  0 

as   n  -»  -   for  all   (x,k)  .   Here   T^  g(x,k)  =  E  (g  (Xfc)  |  XN  (0)  =  x, 


V  (0)  =  k) .   The  presence  of  the  variable   k   prevents  this 
n 


from  being  done  directly.   The  method  we  use  is  to  construct 

co 

a  convenient  sequence  of  functions   <g  >  _,   which  converge 
in  some  sense  to   g.   We  write 

n  oo  ,,      °°  °°    M  M    n  n      ii 

UTjg  -   Ttg|ln    <    ||(Ttg)n   -   Ttg|ln  +    II Tfcg  -   T^ 

+    H^gn   -    <T;g)nlln  (2.6) 


where        refers  to  the  sup  norm  over   S  .   Both   T.   and 
"n  n  t 

oo  , 

T   are  contraction  semigroups. 

oo        oo  , 

<(T^q)  >   n   is  the  sequence  of  functions  constructed 
t^  n  n=l  ^ 

oo 

from  T  g.   Our  goal  is  tb  show  that  each  of  the  three  terms 
on  the  right  side  of  (2.6)  converges  to  0.   The  first  and 
second  terms  can  be  handled  similarly.   For  any  function   g, 

oo 

we  must  quarantee  that  the  constructed   <q  >  -,       sequence 
3  n  n=l     a 

,,  ,.  ..     oo  oo 

satisfies    g   -  g   ■>  0 .   It  will  follow  that    (T.g)  -  T,_g   -►  0 
"^n    ^  n  t^  n    t^  n 

Moreover,  since   ^T4-'  t  _>  0}   is  a  contraction  semi-group 

II T.g  -  T^_g     <  II g  —  g     which  also  converges  to   0.   The 
V3  t^n'n  —  "^    ^n  n  ^ 

oo 

sequence   <g  >  _,   will  be  chosen  in  such  a  way  that  the  third 
term  converges  to  0 . 

We  focus  on  a  convergence  determining  class  of  functions 
g,  those  which  are  bounded  and  have  three  bounded  derivatives. 
For  such  a  function   g(x)   we  define 


gn(x,k)  =  g(x)  +  —  u(x,k)  +  |  v(x,k)     (2.7) 

/n 


where   u   and  v  have  two  bounded  derivatives  in   x   for 

each  fixed  k.   The  functions  u  and  v  will  be  determined 

explicitly  later  and  are  chosen  to  control  the  third  term  in 

(2.6).   Clearly  when   g    is  defined  by  (2.7),  ||g  -g||   =  o(n 

and  therefore  converges  to  0  as  required. 

One  can  apply  the  generator   A    to   g    to  derive 

n      ■Jn 

r 

nQq(x)  +  n1//2[Qu(x,k)  +  g'  (x)  (5  -n  (c+v-k)  )  ] 

+  [Qv(x,k)  +  u  (x,k)  (6  -n  (c+v-k)  ) 

+  |  g"  (x)  (5  -n  (c+v-k)  )  ]  +  0(n"1/2) 


-1/2, 


Angn(x,k) 


if   x  ^  (c+v-k)  //n 

nQg(x)  +n1/2[Qu(x,k)  +  g'(x)(6-  nn1/2x)  ] 

+  [Qv(x,k)  +  uk(x,k)  (6-nn1/2x)  +  \   g"  (x)  (5  +nn1/2x)  ] 


+  Ofn"1/2) 


if   x  =  0,  l//n, ...,  (c+v-k) //n 


where   u  (x,k)  =  -s—  u(x,k) .   Recall  that  Q  annihilates  functions 
of   x   alone,  thus   nQg(x)  E  0.   We  want  to  have   A  g  (x,k) 

converge  to  a  finite  limit  and  to  have  that  limit  be  inde- 

1/2 
pendent  of  k.   For  this  to  occur,  the   n     term  must  be 

controlled  and  the  functions   u   and  v  must  be  chosen  in 

such  a  way  as  to  eliminate  the  variable   k. 


1/2 

The   n     coefficient  in  (2.8)  can  be  rewritten  by 

adding  and  subtracting 

v 

I      tt,  g'  (x)  (5  -n('c+v-k)  )  =  -n  (c+v)  (1-p)  g*  (x)  . 
k=0   K 


We  next  pick   u(x,k)   to  be  a  solution  of 


Qu(x,k)  =-(g'(x)  n (pd-(c+v-k))  +  g' (x)  n (c+v) (1-p)) 
=  -g' (x)  n(k-pv(l-q))  (2.9) 


When   u(x,k)   is  any  solution  of  (2.9),  the  coefficient  of 

1/2 
the   n  '         term  in  (2.8)  becomes 


-g' (x)  n(c+v) (1-p)  if   x  > 


/ 


n 


g'(x)  n((c+v)p  -  n1/2x  -  k)   if   0  <  x  <  c+v"k 

•n 


Equation  (2.9)  can  be  solved  explicitly.   Define 
a,  =  -g' (x)  n(k-P  (l-q))/y,  so  that  (2.9)  can  be  written  as 

K  V 

-Pv(u(x,k)-u(x,k-l))-(k-l) (u(x,k-l)-u(xfk-2))  =  ak_x ,  k  =  1 


.  .  . ,  v 


-v(u(x,v) -u(x, v-1) )  =  a  (2.10) 
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Equation  (2.10)  has  a  solution  since   £i._n  ^v3^  =  0' 
where  ^ir^ir-n      ^s  t^le  stationary  distribution  associated 
with  Q,  or   it,  =  (p „/k*.)/(£._n  P  /il)  .   The  solution  is 
given  by 


u(x,k)  -  u(x,k-l 


k-1 

Jo  "ai     -g'^)^Tk-i 


where 


k 

T   =  I      tt  (1-p  (1-q))   and   T   =  0. 
k    i=Q   IV  v 


Clearly 


k   T 
u(x,k)  =  u(x,0)  -  g'  (x)n  I      ^i^  ,   1  <  k  <  v     (2.11) 

v    .  ,   i-1 


where   u(x,0)   is  arbitrary.   We  let   u(x,0)  =  t  g1 (x)   so 


k    T 

u(x,k)  =  g'  (x)  1  \ 2-  I        -^=i    ,  0  <  k  <  v    (2.12) 

V  2   ypv  i=l   ffi-l  ' 


For  the  choice  of   u   specified  by  (2.12)  we  next  wish 
to  insure  that  the  limiting  generator  is  independent  of  the 
variable   k.   The  function   v   is  chosen  to  eliminate  the 
dependence  on   k.   The  0(1)   term  of  (2.8)  is  given,  for 
x  >  (c+v-k)//n,  by 


11 


Qv(x,k)  +  g"  (x) 


k   T 

I     z1-^   (6  -n  (c+v-k)) 


yp   .",  i.  i 

HKV  1=1    1-1 


+  j  g"(x)  (5  +n  (c+v-k)) 

=  Qv(x,k)  +  H(x,k) • 

Let  H(x)  =  J^_0  TTkH(xfk)   and  consider  Qv(x,k)  +  (H  (x,k) -H  (x)  ) +H  (x) 
We  now  let   v(x,k)   be  any  solution  of 


Qv(x,k)  =  -(H(x,k)  -  H(x))  . 


(2.13) 


Equation  (2.13)  has  a  one-parameter  family  of  solutions, 
since   I^-q  TT,(H(x,k)  -  H(x))  =  0.   When   v(x,k)   is  chosen 
to  be  any  solution  of  (2.13)f  the   0(1)    term  of  (2.3), 
for  x  _>  (c+v-k)//n,  will  become   H(x)   and  will  therefore 
be  independent  of  k.   It  remains  to  calculate   H(x) . 


H(x)  =  g"(x) 


v 
k=0  K 


UP 


k 

I 


Ti-1 


•      1  II   ■      1 

V  1=1     1-1 


(6-n  (c+v-k)  ) 


+  2  (<S  +n  (c+v-k)  ) 


=  cr"  (x) 


5-  r~  I      tt,  (6-n  (c+v-k))  I 


k   T.  , 
l-l 


yp 


v  k=0 


1=1       1-1 


(2.14) 


=  g"  (x) 


5- 


n 


2   v 


V   tt,  (k-p  (1-q))   I 

yp__  ^.Lrs  k        V  L 


k      Ti-1 


v  k=0 


+  n(c+v)(1-p}    y    tt,     y 

yp„   t,£a  k  .£ 


TT  .   , 
1=1     1-1 


k     Ti-1 


k  i=l    ^i-l 
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The  second  term  can  be  rewritten  by  interchanging  the  order 
of  summation.   The  third  term  is   O(l-p) .   We  find 


H(x)  =  g"(x) 


_  _Q_ 
UP 


2  v-1   T. 


v 


X  ^v  L/k^-Pv'1-^' +  o(i-p) 

1=0    i  k=i+l 


=  g"(x) 


2  v-1   T 


6  - 


UP 


I         —  (T  -T.)  +  O(l-p) 


v  i=0 


with   T   =  0   or 
v 


H(x)  =  g"(x)  n 


v-1   T' 


p,  +  — - 
d   pp 


v  i=0 


7T 


+  O(l-p) 


(2.15) 


For  the  functions   u  and   v   specified  by  (2.12 
and  (2.13),  equation  (2.8)  can  be  rewritten  as 


A  g  (x,k) 
n^n 


-n1/2(l-p)  (o+v)  ng*  (x) 


+ 


Pd 


+ 


UP 


v-1   T 

L_  i      ^+  o(i-p) 


v  i=0 


g" (x)  +  o(n  1/2) 


for   x  >  (cfv-k)//n 


=  \ 


(2.16) 


n1/2  n[(c+v)p  -n1/2x-k]g' (x) 


+  ng"  U) 


v-1   T 

Ph  +  HS"  I         ^T   +  O(l-p) 
d    ^pv  i=0   wi 


-  (c+v-k-n    x) 


k  T. 

n    y  _1 

UP    .  r\  1r  ■ 

MKv  i=0  i 


♦  0(n-1/2) 


for   x  <  (c+v-k)//n 
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We  now  introduce  the  "heavy  traffic  approximation." 


In  order  for  the  generator  to  converge  to  a  limiting 

-1/2 
generator  we  must  have   1-p  =  0(n  /  ) .   Specifically,  we 

assume  p  =  p   =  1  -  (9//n)   for  some   6  >^  0  .   In  this  case, 

n1/2(l-p)  =  6,  and  (2.16)  becomes 

r 


A  g  (x,k) 

n^n   ' 


-en (c+v) g1 (x)  +  n 


2-i 

„   v-1  TZn 

^    V    l 

Kd    yp   .  Ln    it. 
v  i=0   i 


g"(x)  +  0(n"1/2) 


for   x  >  (c+v-k)//n 


n[(c+v)p-n  '  x-k]n  /  g'  (x) 


+  ng"  (x) 


2  "i 

n      v-1  TZ 

n   v  i 

p,  +  >  — 

Kd   yp  ■  «  tt. 

v  i=0  i 


-  (c+v-k-n1/2x)  -H-   I    ^ 

ypv  i=0   *i 


for   x  <  (c+v-k)//n 


(2.17) 


We  now  define  a  limiting  generator  A   with  domain 
consisting  of  all  functions   g  having  three  bounded  derivatives 
and   g*  (0)  =  0.   Let 


A  g(x)  =  -6n  (c+v)  g'  (x)  +  n 


P^  + 


2  -i 

v-1  T 

I      — 


yp   •  ft  tt- 
"v  1=0   i 


g" (x) , 


x  >  0 


(2.18) 
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A   is  the  generator  of  a  Markov  process  which 

oo  J  c 

corresponds  to  a  Wiener  process  with  drift   -0n (c+v) ,  scale 


2n 


v-l  T2 

p,  +  — —     )  — 

Kd    yp   >»  it, 

p  v  1=0  1 


-1/2 
and  a  reflecting  barrier  at  0.   The   0(n     )   terms  involve 

the  first  three  derivatives  of   g  which  are  bounded.   It  is 

clear  from  a  direct  comparison  of  (2.16)  and  (2.18)  that 

|A  g  (x,k)  -  A  g|  -*-  0   as   n  -»-  °°   for  all   x  >  0   and   k 
1  n^n         °°^ 

arbitrary.   In  addition,  g*  (0)  =0   is  necessary  for  the 
generator  to  converge  at   x  =  0.   Unfortunately  even  assuming 
g' (0)  =  0, 

2   k   T2 
|A  g(0,k)  -  Aoog(0)|  +  (c+v-k)  -2-  I   ^  g"  (0)   as   n  -*  - 

^pv  i=0   i 

rather  than  to  0.   One  needs  a  special  argument  to  handle 
this  lack  of  convergence  at  the  boundary. 

We  set  out  to  prove  the  third  term  in  (2.6)  converges 
to   0.   Standard  semi-group  results  (see  Burman,  1979,  p.  33) 
give 

<V>„  "  Tt9n  "  I      Tt-S((A~w)n  "  Anwn)dS  (2-19) 

where   w  =  w(t,x)  =  T~g(x) .   Recall  that  wn  =  w  + (l//n) u+(l/n) v 
with   u  and   v   defined  by  (2.12)  and  (2.13)  with  g  replaced 
by  w.    It  follows  that 
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|Tnq    -(T°°q)      II 
1    tyn    v    ty'n"n 


t 
=    11/      t"    -((A  w)    -A  w+A  w-A  w    )dS 


0 


t-S        °°      n      «>        oo        n   n         "n 


£      ^    II^UA^-A^H^dS    +    11^    T^^w-A^dS^ 


t  t 

<       /       ||  (A  w)     -A  w||    dS    +    ||/      t"    e(A  w-A  w    )  dS  || 
—     j.       '      °°     n      °°     n  "^        t-S      oo        n  n        "n 


-1/2 
The  first  term  is  clearly   0(n  '     ) .   It  remains  to 

-1/2 
show  that  the  second  is   0(n    )   as  well.   We  have  shown 

—  1  /2 
|Aw-Aw   =0(n  '     )      except  at  the  boundary  where  it 
oo     n  n 

is   0(1) .   We  split  the  integral  into  two  parts,  for  one 
of  which  the  process  is  away  from  the  boundary,  and  for  the 

other,  near  the  boundary.   The  integral  away  from  the  boundary 

-1/2 
has  an  integrand  which  is   0(n  '    ).   The  integral  near  the 

-1/2 
boundary  is  also   0(n     )   since  under  a  heavy  traffic 

assumption  the  process  is  rarely  near  the  boundary.   The 

details  are  merely  summarized  here;  they  are  based  on  the 

ideas  of  Burman  (19  79) . 

Let   I     be  the  indicator  function  of 
on 


|~o,  ^xds) 


and   I  be    the    indicator   of 

In 
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c+v-k 


We  have 


J   <_S(A„W  -  AnWn)dSHn 


<  I!/'  ^.S(A.W  -  Anwn)IlndS![n  +  11^  T£_S(A.W  -  AnWn)I0ndS||n 


t  t 

<    11/       (A  w   -    A  w    )I,     dS||       +    ||A  w    -    A  w    ||       ||/      T1?    _Irt    dS|| 
—     i,  °°  n   n      In        n  °°  n  n   n      i        t-S    On      "n 


—  1  /2  —1/2. 

The  first  term  is   0(n  /    ) ,  since   |A  w  -  A  w   =  0(n  '    )   off 

00  n   n 

the   boundary.      The    factor      II A  w   -   Aw    ||    =   0(1)  ,    thus    it   remains 
to   show  that 


/      t"   oIn    dS|L    =   0(n    1/2)     . 
~         t-S    On         n 


This    gives    the    total    time    in       [0,t]       spent  near    the   boundary 
We   bound 

by  first  introducing  a  function  h(x)   not  in  the  domain  of 

A^.   We  let  h(x)   have  bounded  support,  be  infinitely 

dif ferentiable  and  be  given  by   h(x)  =  x   for   x  near  0. 

One  can  construct   h  (x)   using  (2.7)  and  apply  A    to 

h    to  find 
n 
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0(1)  if      x    >^  

/IT 

A  h      =     <  (2.20) 

n  n  \ 

n1/2n( (c+v)p-n1/2x-k)    +0(1)       if      x   <   c+v"k 

/n 


One  has 


t 

tV    -    h      =    /      T^A   h    dS 
t  n  n        i        S   n  n 


t  t 

/      t"a  hi.    dS    +    /      T^A  h    In   dS 
^        Snnln  ^        SnnOn 


It   follows    that 


t  t 

/      T^A  h    I„    dS||       <    ||T?h    -   h    ||       +||/      T^A  h    I-    dS|| 
i        S   n  n    On      "n   —  "    t  n        n   n  i         S   n   n   In      "n 


<    2||h    ||       +   0(1)     . 
—      "    n'n 

t 

We   have    shown      ||  /      T^A  h    In    dA||         to   be   bounded   in      n.      An 

11 '         S  n  n  On   "n 

application  of  (2.2)  shows 


||/   t£a  h  In  dA||   =  n1/2nl(C+v)p  -  n1/2x  -  k  +  0(1)  I  ||  /  TcIn  ds  || 
"J    S  n  n  On   "n  K  n 

is  bounded  in   n.   It  follows  that   ||  /   TsIondS"n=  °(n~1//2)- 
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This  finally  concludes  the  argument  which  shows 

i|Ttgn  "  (Ttg)n"n  =  0(n~1/2),  hence  by  (2.6)   ||T*g  -  T~g||n  =  0(n"1/2) 

We  have  thus  shown  that  the  finite-dimensional  distri- 
butions of  the   (X  (t) ,  V  (t))   process  converge  to  those 

n      n  3 

of  a  Wiener  process  with  drift  -0n(c+v)   scale 

v-1  T? 

n  (p  +  J_  I     _L  I, 

d    ypv  i=0  "i 


and  reflection  at  0.   The  diffusion  approximation  treats 
X  (t)   as  though  it  were  such  a  Wiener  process.   For  instance, 
the  limiting  Wiener  process  has  a  stationary  exponential  dis- 
tribution with  parameter 


6(c+v) 


pd    yp  .Ln         i/  1 
v  i=0 


This  is  a  distribution  for   X(nt)//n   and  suggests   X(t) 
will  have  a  steady  state  distribution  given  approximately 
by  an  exponential  with  parameter 


(c+v) (l-p)/(Pj  + 


d   yp 


The  steady  state  mean  data  queue  length  would  then  be 
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2 

n    V_1   T 

P    r    1 

Pj  +  )   — 

d    ypv  i=0  "i 
E(X(t))  =  Y-2lJ1 i_  .  (2.21) 

(c  +  v)  (1-p) 


It  is  interesting  to  consider  the  special  case   c  =  0, 
v  =  1  where  the  two  types  of  traffic  use  the  same  channel. 
Under  heavy  traffic   p  =  pd  +pv/(1+pv),  so   pd  ~  (1+P  ) 
The  mean  data  queue  length  derived  from  the  diffusion  approxi- 
mation (2.21)  will  be 


Pj  +  a_A_W(1_p)^   1  +  a      pv 


"  (i+pv>7/        ~  1_PV     v  d+Pv)2 


The  latter  is  the  exact  expression  derived  by  Fisher  (1978) 
for  this  case.   The  expression  (2.21)  represents  a  generaliza- 
tion of  the  results  of  Gaver  and  Lehoczky  (1979b) .   In  this 
paper,  a  diffusion  approximation  is  given  based  on  the  fluid 
flow  assumption  for  the  data.   For  this  case  the  result  is  the 
same  except  that  the  scale  is  given  by 


2        2 

n  v-l  Tz 


I       — 


rather  than 


yp   •  n     ■""  • 
MKv  1=0    i 


v-l   T2 

d   ypv  i=o  wi 


The  results  derived  in  this  paper  therefore  definitely  generalize 
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the  results  of  Gaver  and  Lehoczky  (1979b)  since  the  variability 
in  the  data  queue  is  now  included.   When  n/y   is  large,  the 
second  term  dominates,  and  the  fluid  flow  approximation  is 
satisfactory. 

The  Wiener  process  approximation  for  the   X(t)   process 
provides  a  method  for  studying  the  dynamics  of  that  process. 
For  instance,  suppose  the  data  queue  were  at  level   x   at 
time  t  where   x   is  large.   One  might  wish  to  study  the  time 
that  elapses  until  the  queue  becomes  empty.   This  is  essentially 
the  duration  of  the  busy  period  under  heavy  traffic  and 
corresponds  to  a  first-passage  time  for  a  Wiener  process.   Let 
us  denote  it  by   T  .   Straightforward  martingale  arguments 
provide  for  its  transform 


-sT 


E(e 


x 


)  =  exp 


'a"'    • 

-  #  - 

\K> 

2 

+ 

2s 

_ 

_J 

(2.22) 


where 


m  =  G(c+v)n  Z     n         (l-p)  (C4v)n 

.2 


n 


,1/2 

n-1  T' 


d    pp 


V 
L 
v  i=0 


It  is  also  easy  to  find  the  mean  first-passage  time 


E  (T  )  =  x/m 
x 


(2.23) 
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One  might  also  be  interested  in  the  area  beneath  the 
sample  path  until  emptiness  occurs,  since  this  area  represents 
the  total  time  waited  by  all  data  customers  involved  in  the 
busy  period.   If   A   represents  this  area,  simple  backward 
equation  arguments  give 


E(V  -  h  +  7?  x  <2-24) 

2m 


2 
where   m  and   a    are  given  in  (2.22) 
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